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ABSTRACT 

We present the analysis of the XMM-Newton EPIC pn spectrum of the Seyfert 2 
galaxy, Mrk 3. We conhrm that the source is dominated by a pure Compton reflection 
component and an iron Ka line, both produced as reflection from a Compton-thick 
torus, likely responsible also for the large column density (1.36^qq4 x 10^^ cm“^) 
which is pierced by the primary powerlaw only at high energies. A low inclination 
angle and an iron underabundance of a factor ~ 0.82, suggested by the amount of 
reflection and the depth of the iron edge, are consistent with the iron Ka line EW 
with respect to the Compton reflection component, being OlOi'igQ eV. Moreover, the 
iron line width, a = 32ji44 eV, if interpreted in terms of Doppler broadening due to 
the Keplerian rotation of the torus, puts an estimate to the inner radius of the latter, 
r = 0 . 6^0 3 sin^ i pc. Finally, two different photoionised reflectors are needed to take 
into account a large number of soft X-ray emission lines from N, O, Ne, Mg, Si, Fe L 
and the Fexxv emission line at 6.71i'iQ'Q2 keV. RGS spectra show that the soft X-ray 
spectrum is dominated by emission lines, while the underlying continuum is best fitted 
by an unabsorbed powerlaw with the same photon index of the primary continuum, 
produced as reflection by a photoionised material with a column density of a few 10 ^^ 
cm“^. We also present the first X-ray spectrum of ROSAT source IXO 30, which 
shows a huge iron line at ^ keV and is well represented either by an absorbed 
powerlaw with F ~ 1.8 or bremsstrahlung emission at a temperature of g keV. 
Its spectral properties point to a likely identification in terms of a weak Galactic 
Cataclysmic Variable, but the lack of any optical counterpart precludes excluding 
other possibilities, like an ULX at the distance of Mrk 3. 

Key words: galaxies: active - galaxies: Seyfert - X-rays: individual: Mrk3 - X-rays: 
individual: 1X030 


1 INTRODUCTION 

Mrk 3 (z=0.0135) is a prototypical Seyfert 2 galaxy, the 
optic al broad lines of which are seen on ly in polarised 
light iMiller fc Goodricblll99ft iTrarJ [l9^ . As soon as it 
was observed in X-rays by GINGA, it was clear that the 
source w as heavily absorbed and presented a strong iron 
Ka line llAwaki et alJllQ^ . |l99l [l. Moreover, a soft excess 
was detected by Ei n stein (lKnrBer_eL-alJ[l99oi) and ROSAT 
iTurner et al ] ll993^ . fIwasawa^^iTi^.994^ suggested that it 
could originate in an extended region, likely dominated by 
scattering of the intrinsic continuum, on the basis of the 
ASGA observation, which also indicated a significant flux 
variability above 4 keV with respect to the previous data. 

* E-mail: Stefano.Bianchi@sciops.esa.int (SB) 


The availability of broad band observations, due to 
RXTE jGeorgmitogoulo^^I 1 11999 D and, in particular, 
B 6000.921 v TcemD^^LFh9T allowed these authors to dis¬ 
entangle three different components: an heavily absorbed 
powerlaw, an unabsorbed Compton reflection component as¬ 
sociated to the iron line and a constant soft X-ray emission. 
The latter component was later found to be spatially ex¬ 
tended along the [O lll] cone and likely produced in a warm 
photoionised material, as shown by the h igh spatial and e n- 
ergy resolution observation by Ghandra iSako et al.ll2nn(i> . 


In this paper, we analyse the XMM-Newton EPIC pn 
and the RGS spectra of Mrk 3, comparing the results with 
the long X-ray history of the source, with the aim to estimate 
some properties of the material around its nucleus. 
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2 DATA REDUTION 

2.1 XNIWI-Newton 

Mrk 3 was observed by XMM-Newton between 2000 Oc- 
tober 19 and 20. w ith all the EPIC C CD cameras, the p n 
iStriider et alJlioOlh and the two MOS iTurner et alJl200ill . 
operating in Full Frame and Medium filter. This paper deals 
only with pn data, the count rate of which is well below the 
maximum for 1 per cent pileup (see Table 3 of the XMM- 
Newton Users’ Handbook). SAS 6.0.0 was used to reduce 
data, with a total net exposure time of 52 ks, after screen¬ 
ing for intervals of flaring particle background, applying the 
pr ocedure to maximise t he signal-to-noise ratio introduced 
bv IPiconcelli et ^ i2004l . An extraction radius of 40 arcsec 
was used for the single and double pattern spectra. After 
a consistency check by fitting the two spectra separately, 
we decided to extract a single spectrum, with all patterns 
from 0 to 4. A weak source is apparent in the pn image just 
outside our extraction region. However, its flux (measured 
from the Chandra spectrum: see below) is completely neg¬ 
ligible with respect to Mrk 3. RGSl and RGS2 spectra of 
Mrk 3 were also extracted, adopting standard procedures, 
with total exposure times of 60 and 58 ks, respectively. 

Finally, another source lies at ~ 1.6 arcmin from 
the nucleus, correspo nding to ROSAT source IXO 30 
iColbert fc Ptakll2002^ . We extracted EPIC pn, MOSl and 
MOS2 spectra for this object, adopting extraction radii of 
19 arcsec, for total exposure times of 47, 49 and 52 ks, re - 
spectively, again after applying the IPiconcelli et alJ (l2004^ 
procedure. In the following fits, the two MOS spectra were 
furthermore added into a single spectrum. Moreover, in or¬ 
der to study the flux variability of this source, we analysed 
all the available EPIC pn observations of the Mrk 3 field 
(listed in Table , with the only exception of observation 
0009220901, where IXO 30 lies partly in the chip gap. 

2.2 Chandra 

Chandra observed Mrk 3 with the ACIS-S HETG between 
2000 March 18 and 19, for ~ 100 ks. The data were already 
analysed bv ISako et alJ J200fil . However, in order to better 
compare them to pn data, we re-extracted first order HEG 
and MEG spectra with ClAO 3.1 and Caldb 2.27, following 
standard procedures. 

The 0th order spectrum of a source ( 0:2000 = 

06'*15™-29A7, (52000 = -k7T’01'40".9) lying at ~ 45 arcsec 
from the nucleus was extracted, to verify if it could con¬ 
taminate the pn spectrum (see above). An extraction region 
with a radius of 3 arcsec was used. The resulting spectrum 
has a 2-10 keV flux of ~ 2.5 x 10“^^ erg cm“^ s“^, thus 
being 200 times weaker than Mrk 3. 

Finally, the 0th order spectrum of IXO 30 
(02000 = 06'*15™-15“, ^2000 = -f71°02'05") was also 

extracted, adopting an extraction region with a radius of 3 
arcsec. 

All spectra were analysed with Xspec 11.3.1. In the 
following, errors correspond to the 90% confidence level for 
one interesting parameter (Ax^ = 2.71), where not other¬ 
wise stated. The cosmological parameters used throughout 
this paper are Hq = 70 km s“^ Mpc“^, Ao = 0.73 and 
go = 0. 


Table 1. Best fit parameters for the continuum (see text for 
details.) 


Nh (cm 2 ) 

1-36/;o p4 X 

r 

1.77 ±0.01 

Ape 

n 09+010 
U.o ^_0 08 

COS i 

0.45“ 

Eo. 6 - 2 keV (erg cm-2 s'l) 

6.3 X 10-13 

E 2 _i 0 keV (erg cm-2 s”!) 

5.9 X 10-12 

Lo.i-isokeV (erg s-l) 

1.1 X 10-14 b 

X^/d.o.f. 

249/208 


“ Xspec default value - ^Absorption-corrected extrapolated lu¬ 
minosity. 

3 DATA ANALYSIS 

3.1 Temporal behaviour 

Figure [i( a) I shows the EPIC pn lightcurve of Mrk 3, along 
with the lightcurve of the background. We can conclude that 
there is no evidence for significant variability during the ob¬ 
servation. Moreover, soft and hard X-rays lightcurves are 
plotted in Figure [I(b)| along with their ratio. Again, taking 
into account sistematic errors due to the subtraction of the 
background, no variability is apparent in the observation. 
As a consequence, we use the whole available exposure time 
in the following spectral analysis. 

3.2 Spectrum 

W e fitted the 0 . 5-13 k eV spectrum with the model adopted 
bvlGapoi et'ai] (Il999ll for the broadband BeppoSAX data, 
which includes a strongly absorbed powerlaw and the com¬ 
ponents arising from reprocessing of this primary contin¬ 
uum from a Compton-thick material, that is a pure Comp- 
ton reflection with the same photon index (model pexrav: 
iMaedziarz fc Zdziarsk^l995^ and an iron Kq line. Moreover, 
a steeper, unabsorbed, powerlaw is used to model the soft 
excess. A column density of 8.46 x 10^° cm“^ was included 
in the model, to take into account the Galactic absorption 
llPickev fc Lockma.n 11199^ ■ 

The resulting fit is unacceptable (x^ = 1074/251 d.o.f.), 
mainly because of large residuals in the soft X-ray spectrum, 
likely due to emission lines. Moreover, the index of the softer 
powerlaw is quite steep (~ 3), so that it is not easy to assess 
its physical origin. Therefore, we examined the RGS spectra, 
before improving the pn fit. We found that the soft X-ray 
spectrum is dominated by emission lines, in particular O vil 
and Ovill Kq lines between 0.5-0.7 keV. Including these 
lines greatly improves the soft X-ray fit of the lower resolu¬ 
tion pn spectrum. We identified each single line with known 
transitions: when the identification was likely affected by a 
blend of lines, we tried to include each line in the model, but 
only when these features improved the reduced ^ pow¬ 
erlaw component is still required by the data, but now its 
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Figure 1. (a) Upper-. Lightcurve of the adopted extraction region for Mrk 3. Middle: Lightcurve for the background, with count rates 
rescaled to the area of the source extraction region. Lower: Background-subtracted lightcurve of Mrk 3. (b) Upper: 0.2 - 2.0 keV 
background subtracted lightcurve for Mrk 3. Middle: 2.0 - 15.0 keV background subtracted lightcurve for Mrk 3. Lower: The ratio of the 
latter to the former. 


photon index can be kept fixed to the one found for the pri¬ 
mary continuum. The final fit is good (x^ = 249/208 d.o.f.): 
the continuum parameters are shown in Table Q while a de¬ 
tailed discussion on the emission line spectrum is deferred 
to Sect. 0 Fig. |2(a^ and p(b)| show the spectrum and the 
adopted best fit model. 

The powerlaw index of the primary continuum and the 
local absorbing column density are well constrained, the two 
values being 1.77±0.01 and I.SGIq q® x cm“^. An inter¬ 
esting result from the fit is an hint for iron underabundance, 
the best fit value being 0. 82~^n'nS (with respect to t he solar 
abundances measured by lAnders fc Grevessel Il98flll . How¬ 
ever, this result must be taken with caution, because this 
parameter is actually calculated on the basis of the iron 
edge depth of the Compton reflection continuum, which de¬ 
pends tightly also on the primary powerlaw index and the 
inclination angle. 

4 DISCUSSION 

4.1 The torus 

The Compton reflection component and the strong iron Ka 
line are clear signatures of reflection from a Compton-thick, 
fairly neutral, material. It is then natural to assume that the 
absorber and the reflector are one and the same material, 
the torus, as usually found in other Compton-thick Seyfert 
galaxies. It is interesting to note that we also detected a 
line consistent with emission from Si less ionised than Si vill 
(see Table 1^. thus likely being another product of reflection 
from the t orus. This featur e was also found in the Chandra 
spectrum llSako et al.ll200fill . 

4.1.1 The Compton reflection 

The 2-10 keV flux of the Compton reflection component only 
is 3.4 X 10“^^ erg cm“^ s“^, completely consistent with the 


3.3 X 10“^^ erg cm“^ s“^ measured by the Chandra obser¬ 
vation, performed seven months before. To check for vari¬ 
ability in past X-ray observations, we can compare the 2-10 
keV total flux reported in Table 0 which by far is domi¬ 
nated by the reflection component (see Fig. |2(b)| . The 1997 
BeppoSA X observation mea sured a flux of 6.5 x 10“^^ erg 
cm“^ s“^ llCaopj et alJIlfifi*^ . slightly larger than the EPIC 
pn one, the difference probably lying in the brighter intrin¬ 
sic co ntinuum of t he sou rce at that time (see below). In¬ 
deed, ICappj et alJ Jl99fl) performed a detailed comparison 
between their dataset and the GINGA and ASCA obser¬ 
vations, concluding that there was evidence of a variable 
component above the iron line energy and a constant one 
below. 


This is also supported by some indications that the 
intrinsic, unabsorbed luminosity (0.1-150 keV) of Mrk 3 
changed between the XMM-Aeruton and Chandra observa¬ 
tions and with respect to the BeppoSAX one, the values 
being 1.1 x 10"^"^, 0.8 x 10"^"^ and 1.3 x 10^^ erg s“^, respec¬ 
tively^. However, it should be stressed that this luminosity is 
driven by the normalization of the strongly absorbed power- 
law, so the measures with X.'MM-Newton and, in particular, 
Chandra are much less reliable than the one obtained with 
the broadband spectrum of BeppoSAX. 


Finally, an estimate of the angle i can be made on the 
basis of the amount of Compton reflection with respect to 
the incident nucle ar continuum which is quite large (i? ~ 1: 
ICappj e~ DlHi), thus requiring an inclination not very 
high to avoid self-obscuration of the torus. 
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(a) (b) 


Figure 2. (a) Best fit model and residuals for the EPIC pn spectrum, (b) The adopted model plotted as E • F{E). 


Table 2. List of the emission lines included in the best fit model, 
originating from a fairly neutral material, likely the torus. 

Energy (keV) Flux“ Id. Et (keV) 


1.744to°J| 0.7±0.1 Sii-viiiKn 1.740-1.769 

6.415 ±0.006 3.8 ±0.2 Fexiii-xvi Ka 6.411-6.424 

7.06* * 0.5 ±0.2 Fei-xviK/3 7.058 


“ 10“® ph cm“^ 

^_Theoret±a]_energies for the transitions (data from NIST: 
l 2 slchenk 2 _g^,gjJ| 222 ^ and references intherein). 

* Fixed 

4.1.2 The iron Ka line 

A prominent iron Ka line is present around 6.4 keV. How¬ 
ever, its EW with respect to the Compton reflection compo¬ 
nent, being GlOitlo eV, falls short of the expected one, which 
should be larger than 1 keV (see e.g. lMatt et alJIlQOGr . Two 
factors can contribute in reconciling the observed value with 
the theoretical one. One is a small inclination angle, since 
the iron line EW increases with this parameter. The other 
is iron underabundance: the dependenc e of EW with thi s 
parameter is almost linear for Afb < 1 llMatt et ;i][l99i. 
Both these requirements are satisfied by the Compton reflec¬ 
tion component properties, which arises in the same mate¬ 
rial: the value of R is consistent with low inclination angles 
(see Sect. 14.1.Ill and the iron edge depth suggests an iron 
underabundance of a factor ~ 0.82 (see Sect. H- 

The iron Ka width is actually resolved in the EPIC pn 
spectrum, with a = 32^^ eV. Even if this measure should 
be taken with care because it is well below the instrument 
resolution, it is interesting to note that it is fully compatible 
with the value found by Chandra, which is ct = eV. 

^ T he latter two value s are those reported by |Sako_ek_alJ ll2000|) 
and|Ca2Ei_£Ls]J (199^, corrected to have the same choice of cos¬ 
mological parameters adopted in this paper. 


The observed iron line width can be in principle explained 
by several means. 

One effect that should be in principle taken into account 
is that the neutral iron line is indeed composed of a doublet, 
Kai a t 6.404 keV an d Ka 2 at 6.391 keV, with a flux ratio 
of 2:1 jBeardenlllOGTll . The difference between the two lines, 
being lower than the spectral resolution of present X-ray 
detectors, is generally i gnored and a w e ighted mean of 6.400 
keV is commonly used T^goob et al.l ll200 Jl already noted 
that a fit with two Gaussian lines does not vary signihcantly 
the measure of the total line width in Chandra gratings data. 
We performed the same test with our high statistics EPIC 
pn spectrum, fitting the neutral iron line with a Kai and a 
Ka 2 , whose energies and relative intensities were freezed to 
the atomic values, while their width was constrained to vary 
together. The resulting, common width is cr = 351^5 eV, in 
full agreement with the value found with a single Gaussian 
line^. 

On the other hand, since the line is produced by repro¬ 
cessing from a Compton-thick material, a Compton Sh oulder 
(CS) is expected on theoretical grounds (lMattll2002ll , as a 
result of Compton scattering of the line photons on the same 
material where they originate. Indeed, the iron line profiles 
of Circinus and NGC 1068 are successfully modelled with a 
narrow core and a CS both in XMM-Aeuiton and Chandra 
spect ra ||Bianclh e^^lj2002 l : lMolendi et all2003l : lMatt et alJ 

I2004l0gle et alJl2003tl . with a better statistical significance 

with respect to a single Gaussian line with a resolved width. 
Therefore, another possibility is that the observed FWHM 
is the resnlt of a bad modelling of the overall profile, which 
does not inclnde the CS. To parametrise this component, we 
added a farther Ganssian line, with centroid energy at 6.3 

^ Note that in order to achieve a fit statistically equivalent to 
the one with a single Gaussian line (whose best fit energy was 
significantly higher than that of neutral iron), we froze the energy 
separation between the two transitions (13 eV, appropriate for 
neutral iron), but allowed the Kai energy to vary. The resulting 
best fit value for the Kai energy is aro und 6.418 keV, i.e. well 
above the one expected for Feix (see e.g. lPalmeri et ^ .1200.31 for 
the theoretical energies of the iron doublet for some Fe ions). 
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keV and cr = 40 eV ||m^^|20o 2)- The CS is not required 
on statistical gronnds (Ax^ — 0 for 1 d.o.f. less), bnt the 
CS flux, even if low, is consistent with production from a 
Compton-thick material, being 10lg% that of the line core, 
and the iron line width is now unresolved (the npper limit 
being 50 eV). We then tested if this model applies well also 
to the Chandra HETG data. However, this is not the case. 
An unresolved line core pins a CS gives a significantly worse 
fit with respect to a single, resolved Gaussian (Ax^ = 13 for 
the same d.o.f.). On the other hand, if the line core width is 
left free to vary (with a best fit value of u = 29lg°), we get 
an upper limit to the CS flux of 13%, still compatible with 
the expectations, but not explaining the FWHM of the line. 

A more likely possibility is that we are not observing a 
line from purely neutral iron. If this is the case, the iron Ka 
line is expected to be actually the blend of a number of emis¬ 
sion features from different ionic species, each composed by a 
doublet, whose transitions are separated by energies of tens 
eV. Indeed, the best fit value of the centroid energy for the 
line in the EPIC pn sp ectrum sugg ests that iron is ionised 
in the range Fe xill-xvi iHous Jl96ai . However, the K(3/Ka 
flux ratio is 0.12^Q'gg. This value does not depend only on 
the ratio between the fluorescence yields of the two transi¬ 
tions (1:8), but also on the different total flux of the contin¬ 
uum and photoabsorption cross sections at the energies of 
the tw o lines. To account for all these effects. iMolendi et alJ 
^2003^ used the lBasl^ il97fjl formulae and obtained for neu¬ 
tral iron a ratio of 0.155-0.160, depending on the inclination 
angle. Therefore, the measured value in Mrk 3, albeit smaller 
than the expected value, is fully consistent with a produc¬ 
tion from neutral iron. Moreover, the K/3 fluorescent yield 
decreases with the ionisation stage, until it becomes null 
for Fexvil, when the M shell is completely void and the 
transition is impossible. Therefore, taking into account the 
above-mentioned value for neutra l iron and scaling it wit h 
the fluorescent yields calculated bv iKaastra fc Mewd lll99,3ll . 
we can estimate that the observed value for the Kfd/Ka flux 
ratio is appropriate for Fex-xi, but is inconsistent with iron 
more ionised than Fexil-xill, just marginally in agreement 
with the centroid energy of the line as observed by the pn. 

Finally, if we assume that the torus rotates around the 
Black Hole (BH) with a Keplerian velocity, it is easy to show 
that the expected FWHM for a line produced in its inner 
walls is approximately 2iik sini ~ 1 300 (Ms/rpc) ^ sini km 
s“^, where the mass is in units of 10® Mq, the inner radius 
in parsec and i is the angle between the torus axis and the 
line of sight. The BH mass of Mrk 3 is estimate d by means of 
stella r velocity dispersion to be 4.5 x 10® Mq iWoo fc IJrrvI 
|20 o 3). This means that the expected FWHM for the iron 
line in this source is 2 770rpc^^^ sini km s“^. Therefore, the 
line width we measure with the EPIC pn, corresponding 
to FWHM=3 hOOj'lj ggg km s“^ if due to Doppler broad¬ 
ening, puts an estimate on the inner radius of the torus: 
r = O.GjjJ'g sin^ i pc. A choice of 30° or 60° leads to central 
values of 0.3 and 0.5 pc respectively for the inner radius. 

Moreover, there is a marginal indication that the line 
flux and energy actually changed between the Chandra and 
the XMM-Newton observation. In Fig. the contour plots 
for the iron line energy vs. flux are shown for both spectra, 
showing that the two curves are not consistent with each 
other at the 90% confidence level (for two interesting pa- 



6.370 6.380 6.390 6.400 6.410 6.420 6.430 

Energy (keV) 

Figure 3. Energy vs. flux contour plot for the ‘neutral’ iron Ka 
line. 

rameters). Even if the reason for this discrepancy is more 
likely to be found in a less than perfect calibration of the 
two instruments (as it may be the case in other high statis¬ 
tics EPIC pn spe ctra of AGN, see e.g. iMolendi et aljE003l : 
iMatt et alJl20oJl . it may be a good exercise to understand 
if such a variation is indeed possible, by estimating the ex¬ 
pected photoionisation time for the neutral iron present in 
the torus. 

The photoionisati on time can be expre ssed in a sim¬ 
plified way, as done bv iRevnolds et ^ (^9^. The formula 
they estimated for O vil can be used also in our case, rescal¬ 
ing for the photoinisation cross section and the threshold 
energ y appropriate for neutral iron fsee lVerner fc Yakovlev! 
Il995ii . so that: 

tph — 1200i?igZ/43 S (1) 

where R = 10^® Rie cm is the distance of the material 
from a source with ionising luminosity {E > 13.6 eV) of 
L = 10"*® Z /43 erg s“®. Adopting an approximate ionising 
luminosity of 10'*'* erg s“* (see TableQl, we find that, for a 
distance of the torus of 1 pc, the resulting tph is of the order 
of 10^ s, i.e. approximately 4 months. This means that a 
variation of the photoionisation structure of the inner wall of 
the torus is perfectly reasonable in a time span of 7 months, 
which separates the two observations. 

4.2 The ionised reflectors 
4-2.1 The soft X-ray spectrum 

The analysis of the RGS spectra shows that the soft X- 
ray EPIG pn spectrum of Mrk 3 is dominated by emission 
lines, in particular from O vil and Ovill Ka between 0.5-0.7 
keV (see Fig. BJ. Therefore, we included these lines in the 
fit with the pn data and added a further powerlaw compo¬ 
nent to model the soft X-ray continuum. There is no need 
now for a very ste ep powerlaw, as was found in BeppoSAX 
iCappi et alJlOo'^ and in our preliminary analysis of the pn 
spectrum. Instead, a powerlaw with the same photon index 
of the primary continuum best fits the continuum underlying 
both the RGS and the EPIC pn emission line spectra. This 
component is most likely due to scattering of the primary 
continuum by a photoionised, Compton-thin gas. Indeed, an 
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Table 3. List of the emission lines originating from an ionised material included in the best fit model, along with their fluxes and their 
likely identifications (see text for details). 


EPIC 

pn 

RGS 




Energy (keV) 

Flux“ 

Energy (keV) 

Flux“ 

Id. 

Et (keV) 

- 

- 

0 43+0-03 

9 7+1.9 
^•'-1.2 

N VI Ka 

0.420 (f) 
0.426 (i) 
0.431 (r) 

- 

- 

b.OU_o 03 

l-4lo,4 

N VII Ka 

0.500 

U.ODO_o 004 

14 4+1-0 

0.562 ±0.001 

0.57 ±0.01 

2 5+1-* * 

Ovii Ka 

Ovii Ka 

0.561 (f) 
0.569 (i) 
0.574 (r) 

U.ODi_o 004 

O-^-0.4 

0.655 ±0.005 

0.666* 

q 9 + 0.6 
'^•^-0.8 

1^-O_0.4 

Oviii Ka 

O VII Kf3 

0.654 

0.666 

n yoc+O.OOS 

U. <'JU_Q QQg 

q '7+0.5 

'^•'-0.3 

0 741+0-005 

9- ''*-L_0.008 

2.8t?l 

O VII RRC 

> 0.739 

0.780t°o°Ti 

1 9+0-4 
^•^-0.3 

0.775* 

0.8t°l 

Oviii K/3 

0.775 

0.822 ± 0.006 

3.2 ±0.3 

0.827 ±0.002 


Fexvii L 

0.826 



0.88+°;°i 

1.6 ±0.4 

O VIII RRC 

> 0.871 

U.»U‘±_Q QQ2 

6.4 ± 0.3 

0-9ltS:o3 

1 S+0-4 
^•^-0.8 

Neix Ka 

0.905 (f) 
0.915 (i) 



0.924 ± 0.002 

1 S+0-4 
^•^-1.0 

Neix Ka 

0.922 (r) 

0.973 ± 0.009 

1 4+0-3 
-*-•^-0.2 

l.UU_Q Q2 

1 4+0.5 
1-^-0.8 

Fexx L 

Fexxi L 

0.965 

1.009 

1.031 ± 0.005 

2.5 ±0.2 

-L-uz_o 01 

1 9+0-4 

Nex Ka 

1.022 

-L-oa_o.o3 

0.4 ± 0.2 

1 n^+O-Oi 
4-1 Jo_0.03 

1 0"*"*^'^ 
l-'^-0.6 

Fexxii L 

1.053 

-‘-•-LU_o.oi 

0 

^•^-0.1 

1.128* 

< 0.5 

Fexxiii L 

1.128 

1.16 ±0.01 

0 9+0-2 

1.170* 

0.8«;« 

Fexxiv L 

1.170 

. o.c+O.OlO 

i.zio_o 007 

1.0 ±0.1 

^■^^-0.04: 

0.8+°« 

Nex K/3 

1.211 

1 310 + 0 - 00 'i" 

1.0XU_Q QQg 

1.3 ±0.1 

1 34+0.01 
J--'3^_0.04 

1-'50_q Q4 

0-8-0.3 

0-+0.3 

Mg XI Ka 

Mg XI Ka 

1.331 (f) 
1.344 (i) 
1.352 (r) 

1.44 ±0.01 

0.7 ± 0.1 

1 47+0.01 
_0.02 

0.8t°;« 

Mg XII Ka 

1.472 

1 r'T+0.02 

-*-■^'-0.03 

0.3 ±0.1 

- 

- 

Mg XI K/3 

1.578 

1 Q/:?f?+0.012 

i.ObD_Q QQ3 

0.6 ±0.1 

- 

- 

Si XIII Ka 

1.840 (f) 
1.855 (i) 
1.865 (r) 

1 99+003 
-^•^^-0.02 

0.2 ±0.1 

- 

- 

Si XIV Ka 

2.005 

9 q'7+0.04 

^■ol _o .03 

0.4 ± 0.1 

- 

- 

Si XIV K/3 

2.376 


“ 10-^ ph cm-2 s-i 

^ Theoretical energies for the transitions. For He-like ions, forbidden, intercombi nation and resonant lin es are separately listed, while 
the photoionisation threshold is listed in the case of RRC (data from NIST: iRalchenko et a ill200l and references intherein) 

* Fixed 
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estimate of the column density of this material can be calcu¬ 
lated on the basis of the ratio between the fluxes of the reflec¬ 
tion component and the primary emission, which is ~ 1%. 
Assuming a covering factor of 0.5, just to have an order of 
magnitude estimate, the column density of the photoionised 
ga s would be approxima tely 3 x 10^^ cm~^. As alread y noted 
bv llwasawa et al and iGuainazzi et al 1 <2004 but see 

also references intherein), it may well be that low resolution 
spectra of highly absorbed Seyfert galaxies are often affected 
by a blending of strong emission lines which mimic a contin¬ 
uum component, and only high resolution spectroscopy can 
unveil the real nature of this kind of soft excess. 

TableOlists all the emission lines likely produced in this 
material which were included in the fit, both in the RGS and 
in the EPIC pn spectra. All the line widths are unresolved, 
so they are assumed as 5 functions. Along with lines iden¬ 
tified with H-like emission from N, O, Ne, Mg and Si, lines 
compatible with blends of He-like complexes of the same 
ionic species are detected. However, the EPIC pn moder¬ 
ate energy resolution does not allow to resolve any of them, 
while only the Ovil, Neix and Mgxi systems are actually 
separated into their forbidden and resonant components in 
the RGS spectra. In these cases (taking tentatively also into 
account the EPIC pn centroid energies of the blends), the 
results suggest an important role of the forbidden line, which 
is usually interpreted as the signature of emission predom- 
inant ly from photoionised ra ther than collisional plasmas 
(e.g. [Porauet fc Dubaull200fll . Moreoyer, the yery presence 
of the O VII and O vill Radiatiye Recombination Continua 
(RRC) argues against a collisional plasma, since in that case 
this feature would be much broader and yer y difficult to de¬ 
tect iLiedahJll99fll : lLiedahl fc PaereljllQQfJl . 

In the Chandra HETG data, some of the He-like com¬ 
plexes are resolved, showing possibly a larger contribution 
from the resonant line (except for O vil) that was interpreted 
in terms of photoexcitation of reso nant transitions, e xpected 
to occur in photoionised plasmas iSako et alJl2nnflll . These 
authors argued against the presence of an additional colli- 
sionally ionized plasma component (that could explain the 
anomalous strong resonance lines) because strong Fe L-shell 
emission is lacking, a result which appears to be confirmed 
by the XMM-Newton spectrum. 

Therefore, plasma diagnostics with both Chandra and 


Table 4. List of the emission lines included in the best fit model, 
originating from a highly ionised material. 


Energy (keV) 

Flux“ 

Id. 

Et (keV) 

6 71 

O- ‘-*-- 0.02 

0.4 ±0.2 

Fexxv Ka 

6.637 (f) 
6.675 (i) 
6.700 (r) 

6 . 97 * 

< 0.1 

Fexxvi Ka 

6.966 

7.60 ±0.05 

0 

^•^-0.2 

Nixxi-xxiv Ka 

7.550-7.668 


“ 10-5 ph cm-2 s-l 

5 Theoretical energies for the transitions. For He-like ions, 
forbidden, intercombinat ion and resonant line s are separately 
listed (data from NIST: |Ralchenko_et_^k||200^ and references 
intherein) 

* Fixed 

XMM-Aeuiton seems to confirm that emission from a pho¬ 
toionised material dominates the soft X-ray spectrum of 
Mrk 3. The (tentative) difference between the role of for¬ 
bidden and resonance lines in the two spectra may be due 
to the spatial resolution of the two instruments, which could 
also play a relevant role, given that the Chandra observation 
showed extended emission which is not spatially resolved by 
XMM-Newton. The XMM-Newton data might simply indi¬ 
cate the importance of photoionisation not only in the in¬ 
nermost nuclear region (as shown by Chandra data) but also 
at larger radii, corresponding to the extended emission re¬ 
vealed by Chandra. In the nucleus, low column densities may 
be the origin of stronger resonant lines, while larger column 
densities in the extended region co uld be responsible fo r pre¬ 
dominant forbidden lines (see e.g. iBianchi et ahllioOSL for a 
full discussion on the relative role of the He-like transitions 
as a function of the column density). 

4-8.2 The ionised iron lines 

A Fexxv Ka line is required in the EPIC pn spectrum, 
with a centroid energy of 6.711 q q 2 ^eV (see Table 3J- This 
line is actually composed of four transitions, the resonance 
line (ui: 6.700 keV), two intercombination lines {x and y: 
mean energy 6.675 keV) and the forbidden line {z: 6.637 
keV). The best fit energy in our data suggests that the dom¬ 
inant transition is the resonance w. This is generally taken 
as a sig n that the gas is mainly in collisional equilibrium 
(see e.g. IPorouet fc Dubaull200(i) . However, the w line can 
be significantly enhanced by resonant scattering, which is 
expected to occur in photoionised plasma. This process is 
very effective at low column densities and the resulting res¬ 
onance line b ecomes the strongest in the He-like ir on spec¬ 
trum (see e.g. lMatt et al.iri996tlBianchi et alJl2005l . and ref¬ 
erences therein). In NGC 1068, the Fexxv w and 2 ; lines 
were actually resolved a nd the large ratio be tween the for¬ 
mer to the latter allowed (Bianchi et alJ ll2005l) to estimate a 
column density of a few 10^^ cm”^. A similar estimate can 
be done for Mrk 3, if the centroid energy of the blend of 
the He-like iron lines is indeed interpreted as the result of a 
dominant w line. 
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Energy (keV) 

Figure 6. The EPIC pn and MOS1+MOS2 spectrum, best fit 
model and residuals for IXO 30. 

A line from highly ionised Ni is also detected, possibly 
associated with the same reflector, while we only found an 
upper limit for a Fexxvi line (see Table A much stronger 
He-like iron line suggests that the ionisation parameter of 
the g as is likely lower than log Ux ~ —0.5 (see lBianchi et all 
I 2 OO 5 I . for details). On the other hand, the ionisation pa¬ 
rameter appropriate for the observed Fe xxv is much higher 
than the one consistent with the production of the other 
lines found in the soft X-ray spectrum, requiring at least 
two different ionised reflecting materials in the circumnu- 
clear region of Mrk 3, similarly t o what found, for example, 
for NGC 1068 iMatt et alJl200'3 . and references therein). 


5 IXO 30 

As shown in figure]^ a bright source is apparent ~ 1.6 ar- 
cmin southeast of the nucleus. The much better astrometry 
of Chandra allows us to locate it at Q 2000 = 06^ 15"“ 15®, 
^2000 = -|-71°02'05", fully con sistent with ROSAT source 
IXO 30 llColbert fc Ptak^^2002^ . 

This source was fir st de tected by ROSAT a nd discussed 
bv lTurner et alJ ^199,3^ and iMorse et alJ {199^. We present 
here the first X-ray spectrum for IXO 30. A combined pn- 
MOS fit with a simple absorbed powerlaw leads to an ac¬ 
ceptable fit = 86/75 d.o.f.), but the addition of an unre¬ 
solved (cr < 600 eV) gaussian line at 6.5 ^q’ 2 keV is required 
at the 97% confidence level, according to F-test. The flux 
of this feature, likely identified as an iron line, even if the 
large errors prevent us from stating whether it is ionised or 
not, is (8 ± 6) X 10“^ ph cm“^ s“^, corresponding to a very 
large EW of 1.0 ±0.8 keV. The final fit is good (x^ = 76/72 
d.o.f.), with Nh = (1.7±0.4) x 10^^ cm“^ and F = 1.77lg;jQ 
(see figurel^. A fit of comparable quality (y^ = 77/73 d.o.f.) 
is achieved by a bremsstrahlung component at a tempera¬ 
ture of 7.5l^ g keV, absorbed only by the Galactic column 
density, plus an iron line with the same properties as for 
the previous model. Finally, we have also tried a single tem¬ 
perature plasma component (model mekal, assuming Solar 
metallicity), which gives a good fit (x^ = 83/76 d.o.f.) with 
a temperature, 8.21^ 2 keV, which is enough to take into 
account the iron line without any further component. 

The observed 2-10 keV flux (for all the models) is 
~ 8.5 X 10“^^ erg cm“^ s“^, corresponding to a luminos- 



Figure 7. Upper: Lightcurve of the adopted extraction region 
for IXO 30. Middle: Lightcurve for the background, with count 
rates rescaled to the area of the source extraction region. Lower: 
Background-subtracted lightcurve of IXO 30. 


Table 5. Log of all the EPIC pn observations of IXO 30 analysed 
in this paper. See text for details. 


Obs. ID 

Date 

Exp. Time (ks) 

1 

2000-10-19 

53 

2 

2001-03-12 

2 

3 

2001-03-20 

10 

4 

2001-03-28 

5 

5 

2001-04-05 

3 

6 

2002-03-10 

5 

7 

2002-03-25 

5 

8 

2002-04-18 

1 

9 

2002-09-16 

5 


ity of around 3.5 x 10'*° erg s“*, if IXO 30 lies at the 
same dist ance of Mrk 3. in g ood agreement with what re¬ 
ported bv lTurner et al.l 11199311 for the ROSAT observation. 
We found a somewhat higher flux, ~ 1.2 x 10“*° erg s“*, 
in the Chandra spectrum of the source, but the fit param¬ 
eters are loosely constrained by the poor statistics, being 
Nh = (1.9/lj'®) X 10°* cm“° and F = 1.9/lo 3 power- 

law model, with an upper limit of 1.3 keV to the EW of the 
iron line. On the other hand, the analysis of all the avail¬ 
able EPIC pn observations does not display any significant 
variability of the source flux within a time span of almost 
two years (see Fig. |^and Table EJ. However, it is interest¬ 
ing to note that the only observation long enough to allow 
some spectral analysis (observation 3 in Table EJ suggests a 
possible variation of the absorbing column density, since a 
larger value, A^h = 7.2^^'^'^ x 10 °* cm is required by the 
data. Finally, no clear evidence for short-term varibility is 
present in the EPIC pn lightcurve of IXO 30 (see figure EJ, 
considering that the high background periods (filtered for 
the extraction of the spectra) are comparable to the source 
flux and are not easily subtracted from the lightcurve. 

We have searched for an optical counterpart of IXO 30 
in the Digitized Sky Survey, but we did not find any possi¬ 
ble companion neither in the POSS-II Red nor in the Blue 
plate, so implying that the source is fainter than R ~ 21 and 
Bj ~ 22.5. It is therefore quite difficult to guess its real na¬ 
ture. A possible interpretation is in terms of an AGN, which 
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Figure 5. Left: EPIC pn image of the Mrk 3 field. Right: Chandra image of the rectangular region shown aside, together with the 
identification of Mrk 3 and IXO 30. 


n-1-1-n 


2 3 4 5 6 7 

Obs. ID 


Figure 8. The history of the observed fiux of IXO 30 among 
the available EPIC pn observations (listed in Table 13. All spec¬ 
tra were fitted with an absorbed powerlaw, adopting the best fit 
values of the longest observation. 


must be a background object, otherwise its luminosity would 
be too low. In any case, the presence of an iron line allows us 
to put an upper limit on the redshift of the source, assum¬ 
ing that the line comes from H-like iron: this limit, being 
approximately « < 0.1, implies that the 2-10 keV luminos¬ 
ity of the source cannot be larger than ~ 2 x lO'^^ erg s“^, 
quite low for an AGN. Moreover, the EW of the iron line, 
though determined with large uncertainty, would be typical 
of a Compton-thick object, but this would be at odds with 
the observed photon index. 

On the other hand, if associated with the Mrk 3 galaxy, 
its luminosity is comparable with that usually found for Ul¬ 
tra Luminous X-ray sources (ULXs) often detected in the 
vicinity of an AGN and p ossibly related to Interm ediate- 
Mass Black Holes (see e.g. iGolbert fc Millerl l20o3l . If the 
latter interpretation is correct, at least a mass of ~ 300 
Mq is required, if the source is emitting near the Edding¬ 
ton limit. However, the origin of the huge iron line remains 


obscure, even if at least another ULX was recently found in 
M82, showing a similar feature iStrohmaver fc Mushotzj^ 
l2003!l . Furthermore, a separation of 1.6 arcmin at the 
redshift of Mrk 3 would imply a large distance of the 
source from the nucleus, being ~ 25 kpc, which is very 
close to the optical 02 ^ diameter of the host galaxy 
ide Vaucouleurs et alJll99lh . so that the probability of con¬ 
tamination f rom background/foreg round objects is signifi¬ 
cant (see e.g. IPtak fc Golbertll2004f) . 

Finally, it is possible an interpretation in terms of a 
Galactic source, like a Cataclysmic Variable (CV). This 
scenario is favoured by the spectral analysis, since a 
bremsstrahlung component and a huge iron line are gen- 
erally the main ingredie nts to model these objects (see e.g. 
iMukai fc Shiokawall993|l . At a distance of ~ 500 pc from the 
Earth, a 2-10 keV luminosity of ~ 2 x 10^° erg s~^ would 
make I XO 30 a weak CV, but s till in the observed range 
(see e.g. lMukai fc Shiokawalli993^ . On the other hand, if we 
assume, as it is generally found, that the secondary of this 
system is a typical M red dwarf, we can get a lower limit on 
the distance of the CV, based on the upper limit we have 
on its apparent R magnitude. The resulting distance, ap¬ 
proximately d > 4 kpc, is in agreement with an origin of 
the source in the Galactic disc, taking also into account its 
Galactic coordinates. This would imply a 2-10 keV luminos¬ 
ity larger than 10®^ erg s“^, which is within the observed 
range in CVs. 


6 CONCLUSIONS 

We have analysed the first XMM-AeM;ton spectrum of 
Mrk 3. The analysis confirm previous results, showing a 
spectrum composed by three main components: a strongly 
absorbed powerlaw, a pure Compton reflection with the 
same photon index associated to an iron Ka line, both pro¬ 
duced as reflection from a Compton-thick torus, and an un¬ 
absorbed powerlaw, again with the same photon index of the 
primary continuum, associated to a large number of emis- 
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sion lines, likely produced as reflection from a Compton-thin, 
photoionised material. 

The iron Ka line EW with respect to the Compton re¬ 
flection component, being only 6 IOI 5 Q eV, is consistent with 
a low inclination angle and an iron underabundance of a fac¬ 
tor ~ 0.82, properties indipendently derived for the torus, 
via the amount of Compton reflection and the depth of the 
iron edge. The iron line width is actually resolved in the 
EPIC pn spectrum, with a = 32j^j4 eV, corresponding to 
FWHM=3 500tl5^g km s"^ if produced by Doppler broad¬ 
ening, thus putting an estimate to the inner radius of the 
torus, r = O.OlJ'g sin^ i pc. Moreover, a possible variation 
on a timescale of ~ 7 months of the ionisation stage of iron 
between the Chandra and the XMM-A^euiton observations is 
fully compatible with the photoionisation time for iron at a 
distance of around 1 pc. 

The soft X-ray spectrum of Mrk 3 is dominated by emis¬ 
sion lines of H- and He-like from the most abundant met¬ 
als, superimposed over an unabsorbed powerlaw with the 
same photon index of the primary continuum. It is impor¬ 
tant to note that the full resolution of these lines in the 
RGS spectra was required before correctly fitting this part 
of the spectrum in the lower resolution pn data, thus pre¬ 
venting the adoption of a much steeper (E ~ 3) powerlaw, 
whose physical interpretation would have been less straight¬ 
forward. From the ratio between the fluxes of the reflected 
component and the primary continuum, a column density of 
a few 10 ^^ cm~^ can be derived for the photoionised mate¬ 
rial. However, it must be noted that at least two different 
ionised reflectors are needed to take into account the emis¬ 
sion lines from lighter metals and the Fexxv emission line 
at keV. 

We have also presented the first X-ray spectrum of 
ROSAT source IXO 30, which shows a huge iron line at 
fl.Slg’j keV and is well represented either by an absorbed 
powerlaw with E ~ 1.8 or bremsstrahlung emission at a 
temperature of keV. Even if the lack of any optical 

counterpart precludes excluding other possibilities, such as 
an ULX at the distance of Mrk 3, its spectral properties 
point to a likely identification in terms of a Galactic CV. 
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